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' Abstract. This work is the sequel of a previous investigation of stationary and cylindrically 

symmetric vortex configurations for simple models representing an incompressible non-relativistic 

Tj ' superconductor in a rigidly rotating background. In the present paper, we carry out our analysis 

with a generalized Ginzburg-Landau description of the superconductor, which provides a prescription 
for the radial profile of the normal density within the vortex. Within this framework, it is shown 

OO ' that the Bogomol'nyi limit condition marking the boundary between type I and type II behavior is 

unaffected by the rotation of the background. 
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. . The present work is the sequel of a previous investigationj on the energy of vortices in a rotating superconductor, 

where it was shown in particular that the magnetic and kinetic contributions to the energy density that are proportional 

to the background angular velocity remarkably cancel. The motivation which brought us to study vortices in rotating 

^ • superconductors is the study of the interior of neutron stars, more specifically their inner core which is believed 

i-rt to contain a proton superconductor. Although we worked out the jpacroscopic description of an array of magnetic 

^ ■ vortices and superfluid vortices in a general relativistic frameworkH necessary for refined analysis of neutron stars, 

O ' thus generalizing the earlier work of Lindblom and MendelEI in a Newtonian approach, we have here restricted our 

I ^1 . analysis to a Newtoniaruframework for simplicity. 

In our previous work,El we left undetermined the explicit structure of the vortex core. One of the purposes of the 
T— I , present work is to provide a specification for the profile of the condensate particle density, based on a Ginzburg-Landau 
K>" ' type approach. Our treatment will however not be restricted just to the widely used standard Ginzburg-Lar 
^^ [ description, but will also be valid for a generalized version (more readily justifiable by heuristic considerations 
leaving arbitrary the coupling constant gc that enters the gradient energy density. _ 

The second purpose of this work is to reexamine the question of the Bogomol'nyi limitQ in the context of a su- 
perconductor in a rotating background. The conclusion will be that the usual boundary between type I and type II 
^\ ' superconductors remains unaffected by the rotation of the background. 

0^ I Before entering the details of this work, let us recall the essential features of the model and define the relevant 

quantities. The superconducting matter will consist of a charged superfluid component, represented by a locally 
2 , variable number density Ug of bosonic particles characterized by an effective mass m, a charge q, and of an ordinary 
r^ ' component of opposite charge which locally compensates the charge of the first component. The essential property 
that distinguishes the superfluid constituent from ordinary matter is that its momentum is directly related to the 
phase variable (/? (a scalar with period 27r) of the boson condensate, according to the expression 



I 

o . 

Q. mv + qA = h\/ ip . (1) 



. i-H , In this formula, h is the Dirac-Planck constant, A is the magnetic vector potential and v the velocity of the Bose 
^^ condensate. The whole system will be described, in addition to the relation (IT]), by the Maxwell equations. 



V X B = 47rj', (2) 

where B is the magnetic field, related to the magnetic potential vector by the usual relation B = W x A, and j* is 
the electric current, which consists of the sum of the currents due to the condensate component and to the ordinary 
component. The ordinary component current will be supposed to be that of a rigidly rotating fluid. To be able to solve 
the coupled system of equations, one needs a prescription concerning the spatial evolution of the condensate particle 
number density. It will be given by an energy minimization principle, using the energy functional corresponding to a 
generalized Ginzburg-Landau approach. 



The plan of the paper will be the following. In section pi we s hall use the cylindrical symmetry and introduce new 
variables to simplify the system of coupled equations. Section |lll| will be devoted to the energy minimization principle, 
which will give a prescription for the determination of n^ . And finally, section |^ will deal with the Bogomol'nyi limit 
condition. 

II. SYSTEM OF COUPLED EQUATIONS 

The scenarios we shall consider will be of the usual kind, in which each individual vortex is treated as a stationary, 
cylindrically symmetric configuration consisting of a rigidly rotating background medium with uniform angular velocity 
f2oo, say, together with a charged superfluid constituent in a state of differential rotation with a velocity v, which 
tends at large distance towards the rigid rotation value given by ^caf, where r is the cylindrical radial distance from 
the axis. It will be supposed that the superfluid particle number density Us is a monotonically increasing function 
of the cylindrical radius variable r, tending asymptotically to a constant value, rioo, say, at large distances from the 
axis. It will be supposed that the local charge density is canceled by the background so that there is no electric field, 
but that there is a magnetic induction field with magnitude B and direction parallel to the axis, whose source is the 
axially oriented electromagnetic current whose magnitude j will be given by 

j ^ qns{v ~ noor) . (3) 

The relevant Maxwellian source equation for the magnetic field (0) will have the familiar form 

- = ~4nj . (4) 

The other relevant Maxwellian equation is the one governing the axial component A (which in an appropriate gauge 
will be the only component) of the electromagnetic potential covector, which will be related to the magnetic induction 
by 

'-^^rB. (5) 

dr 

The essential property distinguishing the superconducting case from its "normal" analogue is the London flux quan- 
tization condition (j^, which in the present context (where all physically relevant quantities depend only on the 
cylindrical radius r) will be expressible in the well known formQ 

mv + qA — , (6) 

r 

where N is the phase winding number, which must be an integer. 

Before proceeding, it will be useful to take advantage of the possibility of transforming the preceding system of 
equations to a form that is not just linear but also homogeneous, by replacing the variables v, B, A by corresponding 
variables V, B, A that are defined by 

V = v-n^r, (7) 

B = B-B^, (8) 

A^A- irSoo , (9) 

where Bqo is the uniform background magnetic field value that would be generated by a rigidly rotating supercon- 
ductor, which is given by the London formula 

9m 

B^^-—n^, (10) 
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obtained by combining (ph and (m in the special case of rigid corotation, i.e., with v — floor. 
In terms of these new variables the equation (ph will be transformed into the form 



^^.S, (11) 

dr 

while the other differential equation (0) will be transformed into the form 

^^-4,i, (12) 

in which, rewriting equation (0), we shall have 

i = qUsV . (13) 

Finally, the flux quantization condition (|6|) will be converted into the form 

Nh 

mV + qA= , (14) 

r 

which can be used to transform (0) into 

m d{Vr) 



qr dr 



-B. (15) 



The advantage of this reformulation is that unlike w, B, and A^ the new variables V, B and A are subject just to 
homogeneous boundary conditions, which are simply that they all tend to zero as r — + oo. 

III. ENERGY MINIMIZATION PRINCIPLE 

The equations of the previous section are not sufficient by themselves to fully determine the system. In order 
to specify the radial distribution of the condensate particle number density Us we will use an energy minimization 
principle based on a model in which the condensate energy density is postulated to be given as the sum of a gradient 
contribution and a potential energy contribution by an expression of the form 

•^con — ^grad + V , (16) 

where the contribution £grad is proportional to the square of the gradient of Us with a coefficient that, like the potential 
energy contribution V , is given as an algebraic function of Ug by some appropriate ansatz. The use of such a model as 
a fairly plausible approximation is justifiable by heuristic considerationsQ that motivate the use of an ansatz of what 
we shall refer to as the Ginzburg type, according to which the energy contribution is postulated to have the form 

^glh^fdnsV 

igrad - g -T- , (17) 

8m Us \ dr J 

where gc is a dimensionless coupling constant, while the potential energy density V is given in terms of some constant 
proportionality factor £c, say, by the formula 

V^8jl^^)\ (18) 

which provides a particularly convenient ansatz for interpolation in the theoretically intractable intermediate region 
between the comparatively well understood end points of the allowed range < ris < Uoo- The constant £c is 
interpretable as representing the maximum condensation energy density. Its value is commonly expressed in terms of 
the corresponding critical value He, say, representing the strength of the maximum magnetic field that can be expelled 
from the superconductor by the Meissner effect, to which it is evidently related by the formula 

8.^^. (19) 

The total energy density associated with a vortex will be of the general form 



P — z. °° foil 



where i^mag, ^kin and £con are respectively the magnetic, kinetic and condensate energy contributions. More precisely, 
i^mag is the extra magnetic energy density arising from a non-zero value of the phase winding number N , i.e., the local 
deviation from the magnetic energy density due just to the uniform field B^o (associated with of the state of rigid 
corotation characterized by the given background angular velocity fioo), namely, 

while £kin is the corresponding deviation of the kinetic energy from that of the state of rigid corotation characterized 
by the given background angular velocity fJoo, namely, 

fkin = ^ns (w^ - 17oo V2) . (22) 

It is convenient for many purposes to express such a model in terms of a dimensionless amplitude ^ that varies in the 
range < '0 ^ 1 according to the conventional specification 

Us = -fA^^-oo • (23) 

Within the general category of Ginzburg type models as thus described, the special case of the standard kind of 
Ginzburg-Landau model is characterized more specifically by the postulate that the gradient coupling constant gc 
should be exactly equal to unity. This ansatz has the attractive feature of allowing the theory to be neatly reformulated 
in terms of a complex variable ^ = ij}^"^^ where ip is the phase that appears in (Q) in a manner that is evocative of the 
Schroedinger model for a single particle. Indeed, it is easy to verify that in the case of ^c = 1, the gradient term ( p^ ) 
and the kinetic term in (p3) can be rewritten, using dl^) and (p3|), to give the usual Ginzburg-Landau type gradient 
term, i.e., 
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(24) 



where the covariant derivative is defined as 

(25) 

However, although there are physical reasonsQ for expecting that gc should be compasAble with unity, the seductive 
supposition that it should exactly satisfy the Landau condition ^c = 1 is more dubious. cl This more specialized ansatz 
will not be needed for the work that follows, which applies to the generalized Ginzburg category with no restriction 
on the parameter Re- 
using the expression 

Stt 2q Sirr dr 

for the first two terms in the combination (EO) , it can be seen that the total energy density arising from the presence 
of the vortex will be given by 



f mag + f kin = — + —jV + ^ — (r^S) (26) 



6^ i^ 2)2 ^ioo 

Stt Stt ^ ^ '' 2™ 



9cTt---j- ] + {mVipf 



Boo d(r^B) 

-^^-^. (27) 

OTrr dr 



The equation governing the distribution of the condensate particle number density Us is obtained by requiring that 
the integral of the energy density (|2^) be stationary with respect to variation of Ug or equivalcntly of ip, which gives 
the field equation for the latter in the form 

'-^^(r'^)-rnV^i^-—SMl-i^')- (28) 

mr dr \ dr J rioo 



When this equation is satisfied, it can be seen that the energy density (27) will reduce to a value given by 



B^ Hc^ , 1 d f r^B g^n^i^d^P , 



Stt Stt r dr \ 8n Am dr 

in which the last term is a divergence that goes out when integrated, so that for the total energy per unit length one 
is left simply with 

U^^J{l3' + Hc^l~^/))dS. (30) 



IV. BOGOMOL'NYI INEQUALITY 

We shall now try to rewrite the energy density associated with the vortex in a different form. Let us begin by 
writing the relation 

(5en^)V(mV^)^= (,.?.^T-V^)'±..?.#^S±^^, (31) 

which can be obtained by rewriting 'ip{dij)/dr) in terms of dj / dr and dV/dr, the latter term being transformed by use 
of the Maxwell equation (fsh. In the first term on the right hand side of (|3l|), one takes the minus sign if V is positive, 
the plus sign otherwise. We are thus able to obtain a Bogomol'nyi type reformulation of (E^) that is given by 



which generalizes the original versionfl by inclusion of the term with the coefficient B^ , which allows for the effect of 
the background rotation velocity f2oo- Since this extra term is just the divergence of a quantity that vanishes both 
on the axis and in the large distance limit it gives no contribution to the corresponding integral expression, which 
therefore has the same form as the usual Bogomol'nyi relation for the non rotating cascB More generally, performing 
the BogomoPnyi trick {pW just for a fraction / of the combined kinetic and gradient contribution, one can see that 
for a vortex with (relative) magnetic flux 



$= / BdS, (33) 

the total vortex energy per unit length and per unit of flux will be expressible in the form 

u^H^m + ^ii-l^) [{i-^^fds+^jm-l^^^ds 



(,,6^_„.|V|,)'<,S+ <i^/((,.4)%(,nV*)^) iS. m 



in which k is a dimensionless constant given by the deflnition 



2gcTTfiqnoo 
It can be convenient to introduce the so-called London penetration length A, deflned by the expression 



(35) 



and the usual flux quantum 



iTTq^Uoo 



2TTh , , 

$ = . (37) 
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It is to be observed that all the terms in this expression will be non-negative provided the quantity / is chosen not 
only so as to lie in the range < / < 1 but also so as to satisfy / < k, a requirement that will be more restrictive if 
K < 1. In the latter case we can maximize the first term on the right of ( |34|) by choosing f — k, thereby incidentally 
eliminating the second term , so that we obtain the lower limit 

<& 47r 



If K > 1 we shall be able to choose / = 1 , thereby eliminating the final term in (|34|) . so that we obtain the inequality 



«>1 ^ lL>?^- (39) 

More particularly it can be seen that choosing / = 1 will eliminate both the second term and the last term on the 
right of (^) in the special Bogomol'nyi limit case characterized the condition 

K = \. (40) 

(Readers should be warned that much of the relevant literaturdj'u follows a rather awkward tradition in which the 
symbol k is used for what in the present notation scheme would be denoted by ^f2n^ which instead of (^0|) makes the 
critical condition come out to be k = l/\/2-) In this Bogomornyi limit, the energy per unit flux per unit length (p3) 
will be minimized by imposing the conditions 

5,;i^ = m|V|^, (41) 

dr 

(with the sign adjusted so as to make the right hand side positive), and 

which (as when the background is non-rotatingcO) will automatically guarantee the solution of the field equations in 
this case, annihilating the last two terms in (|^ so that one is left simply with 

The qualitative distinction between what are knownliil as Pippard type or type I superconductors on one hand and 
as London type or type II superconductors on the other hand is based on the criterion of whether or not, for a given 
total flux, the energy will be minimized by gathering the flux in a small number of vortices, each with large winding 
number 'N or by separating the flux in a large number of vortices, each just with unit winding number. Within the 
framework of our analysis, a model may be characterized as type I if the derivative with respect to \^\ of C//|*I'| is 
always negative, i.e., if cff//(i|$| < C//|<&|, in which case the vortices will effectively be mutually attractive, and as 
being of type II if the derivative with respect to \^\ of t//|<I>| is always positive, i.e., if df//(i|$| > C//|$|, in which 
case the vortices will effectively be mutually repulsive. One can of course envisage the possibility of models that are 
intermediate in the sense of having a sign for the derivative of f7/|$| that depends on A'^, so that the minimum is 
obtained for some large but finite value of the winding number. 

What can be seen directly from (E^) is that within the category of models characterized by the Ginzburg Landau 
ansatz, the special Bogomol'nyi limit case lies precisely on the boundary between type I and type II, since it evidently 
satisfies the exact equality 

for all values of |$|. The implication of our work is that the well known conclusiorB^Q that there is neither repulsion 
nor attraction between Ginzburg model vortices in the Bogomol'nyi limit case will remain valid even in the presence 
of a rotating background. (It haa,also been shown to be generalizable to cases where self gravitation is allowed for in 
a general relativistic frameworkE3E3) . 

Since the change of variables performed in section O has transformed the system to a representation that is formally 
identical to that of the case with a non-rotating background, it can also be concluded that the conclusions of the 
pioneering analysis of KramerD will remain valid, i.e., that in the sense of the preceding paragraph the system will be 
of the Pippard kind (type I), if k < 1, i.e., 

47rA^7?, < .9,1. ^ |^<^' (45) 

and that it will be of the London kind (type II), if k > 1, i.e., if 

4.A^i7. > .9,$ ^ ^>^' (46) 

6 



at least so long as the winding number N and the parameters k and gc are not too far from the neighborhood of unity, 
in which the numerical investigations have been carried out. It does not seem that the possibility of an intermediate 
scenario, with C//|$| minimized by a winding number N that is finite but larger that one, can occur within the 
framework of Ginzburg type models except perhaps for parameter values that are too extreme to be of likely physical 
relevance. 
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